%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
% Fonction d'assemblage de la matrice de raideur par patch
%
% Entr�es :
%
% sous-maillage : vecteur de coordonnées des noeuds du sous-maillage
% uh : [vecteur 4x1] extraction de la solution éléments finis
% epsilon : =1 sur le premier �l�ment du patch,
%           =-1 sur le deuxi�me �l�ment du patch,
%
% Sorties :
%
% ie : [matrice 4x4] matrice ie

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

function [ K ] = Assemblage_K_patch_fluxfree( Ke1,Ke2 )

n1=size(Ke1,1);
n2=size(Ke2,1);

K=zeros(n1+n2-2,n1+n2-2) ;

K(1:n1,1:n1)=Ke1;
K(n1-1:n1+n2-2,n1-1:n1+n2-2)=K(n1-1:n1+n2-2,n1-1:n1+n2-2)+Ke2;

end

